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We prove that quantum many-body systems on a one-dimensional lattice lo- 
cally relax to Gaussian states under non-equilibrium dynamics generated by a 
bosonic quadratic Hamiltonian. This is true for a large class of initial states — pure 
or mixed — which have to satisfy merely weak conditions concerning the decay 
of correlations. The considered setting is a proven instance of a situation where 
dynamically evolving closed quantum systems locally appear as if they had truly 
relaxed, to maximum entropy states for fixed second moments. This furthers the 
understanding of relaxation in suddenly quenched quantum many-body systems. 
The proof features a non-commutative central limit theorem for non-i.i.d. random 
variables, showing convergence to Gaussian characteristic functions, giving rise 
to trace-norm closeness. We briefly relate our findings to ideas of typicality and 
concentration of measure. 
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1 Introduction 

In what sense can closed local many-body systems in a non-equilibrium situation relax 
to an apparent equilibrium? Instances of that question have a long tradition in the litera- 
ture. This apparent contradiction of having entropy preserved in any closed system and 
at the same time arrive at an equilibrated situation for long times can yet be resolved by 
acknowledging that merely some observables, and in particular entire subsystems, can 
well appear as if they had truly relaxed, even in closed systems. In local observations, 
such a composite system would then look entirely equilibrated. 

Early work formulating a quite similar intuition, albeit not in quantum lattice sys- 
tems, was concerned with fermions freely moving in space [11. The seminal work in 
Ref. |2][3) rigorously developed such an intuition for classical harmonic systems on cu- 
bic lattices. Refs. 0[5]|6) consider that question in instances of a bosonic or fermionic 
fully Gaussian setting, where initial states are Gaussian, and the evolution is governed 
by a quadratic Hamiltonian. Ref. Q then arrives at a true local relaxation theorem for 
subsystems where initial states are not taken to be Gaussian: Yet, in the course of the 
non-equilibrium dynamics of the system, locally, the states become Gaussian, so max- 
imum entropy states for given second moments. This is true without having to invoke 
a time average. In such lattice models, the intuition is that non-equilibrium generates 
local excitations at each site that then travel ballistically through the lattice, resulting in 
a mixing at each site with excitations from further and further separated sites coming 
in HUE). 
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Some of the revived significant recent interest in equilibration and relaxation in 
closed quantum many-body systems, needless to say, has been sparked off by novel 
experiments with cold atoms in optical lattices or on low-dimensional structures. For 
theory work relating to those experiments, see, e.g., Refs. ITT1 [121 [131 [T4l [TUl |9l and 
references therein. This development is quite intriguing, as questions of an apparent 
local relaxation in non-equilibrium dynamics can be measured and probed under very 
precisely tunable laboratory conditions. 

In this work, we present a fully rigorous relaxation theorem, generalizing the find- 
ings of Ref. 0, for one-dimensional quantum lattice systems. Physically speaking, 
the results presented here can be viewed as a relaxation theorem forming a "theoreti- 
cal laboratory", describing an idealized situation of non-equilibrium dynamics in the 
Bose-Hubbard setting: The system being initially held in some ground or thermal state 
of the full Bose-Hubbard-Hamiltonian, and then suddenly switched or "quenched" to 
a parameter regime of strong hopping. 

It is shown that for general pure or mixed initial states satisfying quite weak as- 
sumptions concerning the clustering of correlations — properties that are expected to 
hold in particular for ground states of local many-body models — , locally, the system 
will relax to a Gaussian state under the dynamics generated by a quadratic Hamilto- 
nian. That is to say, we prove the "local relaxation conjecture" for one-dimensional 
clustering states under quadratic Hamiltonians on a ring of L sites. To wit, for an ini- 
tial state go that is evolved in time under a quadratic Hamiltonian and a given finite 
subset S of the chain, we show that for every e > and sufficiently large L there is 
a relaxation time i re iax, and a recurrence time t KC (which grows unboundedly with L), 
such that the reduction of the time-evolved state to S, Qs(t), fulfils 

\\es(t) - QG(t)\\tr < e for all i re i ax < t < t rec , (1) 

in particular, 

lim lim Us(t)-QG(t)\\ b = 0. (2) 

t— >oo L—too 

The state ga{t) is a Gaussian state with the same second moments as Qs(t). The proof 
features the explicit forms of t re iax and t lec in terms of the coupling parameters in the 
Hamiltonian and the decay behaviour of correlations of the initial state. 

The key ingredients in the proof are a quantum Lindeberg central limit theorem, 
a notion of locality of dynamics as is manifest by a Lieb-Robinson bound, and a 
Bernstein-Spohn-blocking argument reminiscent of the classical situation. Essentially, 
we will show that for a given time excitations arising from far away sites will hardly 
influence the dynamics at a subset of sites, whereas the close sites give rise to a mixing 
such that a Lindeberg central limit theorem can be invoked. In several ways, our quan- 
tum treatment reminds of the classical situation presented in Ref. |2 , 3 ], with the block- 
ing argument being essentially identical, with significant differences in many other as- 
pects. In this work, we focus on the important case of a one-dimensional setting and 
allow for assumptions on the initial state that can be relaxed, in order to render the 
argument simpler than a fully general argument lfl5l . 

After rigorously proving convergence of characteristic functions to Gaussian func- 
tions in phase space, we prove closeness in trace-norm of the respective reduced states. 
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In this way, we derive in the trace-norm topology a closeness to maximum entropy 
states under the second moments. In some ways, the arguments reminds of the situa- 
tion of discrete time in convergence to stationary states in quantum cellular automata 
|fl6l . We also finally briefly discuss the findings in the light of recent results concerning 
concentration of measure and notions of typicality ifTTl [THl |T9l . 

2 Preliminaries 

In this section, we are collecting preliminaries that are being used in the proof. We 
also specify the model under consideration here and will specify the assumptions on 
the initial states in the subsequent section. Note that the assumptions on such initial 
states are very mild and do by no means only include pure or Gaussian states; merely 
the model Hamiltonian as such is taken to be a specific nearest-neighbor Hamiltonian. 

2.1 System and dynamical setting 

We consider a chain C = {1, . . . ,L} equipped with periodic boundary conditions. 
Each lattice site is associated with a bosonic canonical degree of freedom equipped 
with the usual symplectic form. The distance between lattice sites i,j € C is due to 
the periodic boundary conditions given by 

dij = mm{\i - j\,\i - j + L\,\i - j - L\} . (3) 

For future purposes, we also define distances between subsets of sites. For A,Bc£ 
we define their distance as 

dA,B = min d i; j. (4) 

Starting from an initial quantum state (neither restricted to be a pure state nor a 
Gaussian state) qq on £, we are concerned with its time evolution, 

g(t)=e- i " t Q e i6t , (5) 

under the Hamiltonian 

A = l E (HAi/bj + hAi/bf). (6) 

bi denotes a bosonic operator at site i. We take the Hamiltonian to be defined by 

A itj = -JS diJtl . (7) 

We will work in units in which J = 1, H = 1. This is the Hamiltonian of a har- 
monic chain, in particular approximating the Hamiltonian of a Bose-Hubbard model 
in the limit of a hopping that is dominant compared to the interaction. Note that lo- 
cal relaxation may be shown for general quadratic Hamiltonians with sufficiently local 
couplings (exponentially decaying in the coupling distance) on (not necessarily cubic) 
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lattices of arbitrary spatial dimension lfl5l . The assumption that H is quadratic, is es- 
sential, however, for the validity of the argument. Yet, the physical intuition that local 
relaxation is due to a mixing of excitations that essentially travel ballistically through 
the lattice should be expected to be valid more generally. In this sense, the present 
proof constitutes the mentioned "theoretical laboratory" in which this mechanism can 
be studied in great clarity. 

In particular, we will consider local properties of g(t), and will see that in a sense 
this reduced state will eventually relax in a sense yet to be specified. To this end let 
S C C be a subset of sites of the lattice and define the reduced state associated with 
the degrees of freedom of this subset as 

&(*)=tr A5 [g(t)]. (8) 

We write its characteristic function x§ s (t) '■ <C as the expectation value of the 

Weyl operator, i.e., as 

X Ss( t)(P)=trs[Qs(t)D(l3)}, D((3) = l[e^-^ : (9) 

as a complex-valued phase space function uniquely defining the state. States the char- 
acteristic function of which are a Gaussian in phase space are referred to as being 
quasi-free or Gaussian states. 

One finds, after solving Heisenberg's equation of motion and using the cyclic in- 
variance of the trace, for the time evolution of the characteristic function of the reduced 
state Qs(t) 

X6s (t){P) =Kc[Qv\[e a ^ m - aUmk ] =:txc[Q D(cx(t,l3))}, (10) 

iec 

where a G (D' £ ', 

= c = e~ itA , (ii) 

the latter to be understood as a matrix exponent of the matrix A collecting the coupling 
coefficients. Note that this is exactly the setting of a sudden quench: One considers a 
ground or thermal state of some Hamiltonian (or any initial state) and the abrupt switch 
to the Hamiltonian H, and follows the free time evolution under this new Hamiltonian. 
Note that the weak assumptions on the decay of correlations also allows for initial 
ground states of critical models. 

2.2 Moments 

For operators O defined on the Hilbert space associated with the lattice model, we write 
expectation values as 

(6)=tx[g d]. (12) 

For vectors of phase space coordinates a, a' € (D' £ we define the anti-hermitian 
operator 

b(a) = J2 («$ - a$U) , (13) 
iec 
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and instances of second and fourth moments of b(a) as 

<r(a)=a(a,a), *(a, a') = (6(0)6(0')), f(a) = ([b(a)] i ). (14) 
For A C C we define the vector a A e (D L by 

(a A )i = i . (15) 

I otherwise, 

as the indicator function and write 

VA,B = v(<XA,Oib), <TA = (?A,A, fA = f{oLA), (16) 
for certain second moments and functions of indicator functions. 

3 Assumptions on initial state and locality of dynamics 

In this section, we state the assumptions on the considered initial states and highlight 
the role of Lieb-Robinson bounds in the argument. 

3.1 Assumptions on the initial state 

The assumptions on the initial state essentially express some degree of clustering or 
decay of correlations. This clustering can be assumed to be quite weak, and even a 
slow algebraic decay of correlations is allowed for. We will also include some natural 
conditions on the initial state that significantly simplify the argument. In particular, this 
includes an assumption of the initial state to commute with the total particle number 
operator. In any bosonic system involving massive particles, this natural requirement 
will always be satisfied. 

More specifically, we assume that the initial state is such tha{] 



0. (17) 



iec 

Then, for all a E <D L and all A,B C Cwe have (b(a)) = 0, and 

°A,B = p{<xa, «s) = - ^2 ^2 aiOL *3 + (Mi)) 1 

ieAjeB 

fA = f(a A ) = ((Wjhk) + (^MYk) + (blhbify (18) 

i,j,k,leA 

+ (SfcStSJg,) + <S fc StSiS|) + (SfcSiSJSt)). 

The first assumption concerns two-point correlations. 

1 This simplifies the calculations significantly as expectation values of products involving unequal num- 
bers of annihilation and creation operators vanish, is, however, not necessary to show local relaxation 1 1 51 - 
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Assumption 1 (Two-point correlations) Let the initial state qq be such that there ex- 
ist absolute constants c\ , \i\ , e\ > such that 



< 



Cl 



[1 + d t J 2+t " 1+C1 



(19) 



for all i, j. 



This is a natural assumption on the initial state in a way that does not depend on the 
system size. This assumption may be relaxed to decay stronger than the spatial di- 
mension of the lattice (i.e., in the setting at hand the two may be replaced by a one), 
resulting, however, in a much more involved proof UTI . We note a first immediate con- 
sequence of Assumption Q] that will be used in the main theorem. We also encounter 
the recurrence time: the system has to be sufficiently large for the bounds on the second 
moments to hold. 

Lemma 1 (Variances) Under Assumption\J}and for L > \t\ 7 / 6 > 1 one has the fol- 
lowing bounds^ 

a mm{\A\, |B|} 

Wa,B\ < C2\\p\\l 



for all A, B C C, 



i=l 



|i| 2 / 3 [l + d A , B ] 1+tl1 



< Wc\a\ + 2Wc\a,a\ + c2\\[3\\IJ2 



\Ai\ 



i=l 



\t\ 2/3 [l+d A \ Ai , Ai ] 



for all Ai C C, i = 1, . . . ,n, A = UiAi, and 

Wa,b\ < c 2 ||/3||22-^ s 
for all A,B C C such that 4e|t| < d^g. Here, 

c 2 = 37 2 (4ci+2)C(l + ei) 
and Ci, ei, and fix are as in Assumption^ 
The second assumption concerns four-point correlations. 



Assumption 2 (Four-point correlations) Let the initial state p be such that there 
exist absolute constants C3 , e 2 > such that 



(20) 



(21) 



(22) 



(23) 



VM k k 



< 



E 



c 3 



([l + dr,.][l + dt,u]) 



1+C2 



(r,s,t,u)EP(i,j,k,l) 

for all i, j, k, I. Here, P{i,j, k,l) contains all permutations of the indices k, I. 



(24) 



2 We use the following vector norms for /3 6 C' s ': ||/3||i = 5TJ;gs 
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Lemma 2 (Fourth moments) Under Assumptions \l\2\ and for L > \t\ 7 / 6 > 1, one 
has 

n 

Ei/^i^^iiipii? 2 ^ 1 (25) 
i=i 

for all Ai C C, i = 1, . . . ,n, vv/f/z .Aj n A, = 0/or i 7^ j. Here, 

c 4 = 96(6c 3 + 3(4 Cl + 1))37 2 C 2 (1 + min{ £l , e 2 }), (26) 
C\, ei are as in Assumption\l\and C3, £2 a* in Assumption^ 
We finally put a clustering condition on the initial state. 

Assumption 3 (Clustering) Let the initial state qq be such that there exist absolute 
constants €5,(12 > such that 



< _5 (27) 

" [i + <Ufl] 1/2+Ma 



;/3(a^)/3(a B )) - (D(a A ))(D(oc B )) 
for all A,Bc£ and all C. 

Note that this condition may again be relaxed lfl5l . 

3.2 Lieb-Robinson bounds and locality of dynamics 

Lieb-Robinson bounds are upper bounds to group velocities of excitations travelling 
through a quantum lattice system. They define a causal cone of a local excitation, out- 
side of which any influence of this excitation is exponentially suppressed. As such, they 
provide an upper bound to the speed of any non-negligible information propagation by 
time evolution of a quantum lattice model l20l |2T1 l22l l23ll . Originally formulated for 
spin systems, in the present context, we need to invoke an instance that is suitable 
for the harmonic, infinite-dimensional individual constituents at hand E4l l25l . Sub- 
sequently we will formulate a useful related bound defining a causal cone for matrix 
entries of the propagator itself, which will be used in the main theorem. 

Lemma 3 (Lieb-Robinson bound) For all times, all i,j e C and all A C C with 

4e|i| < d,A,s one nas 

icy < LL H ^ m\i^ dA - s - (28) 

For all t, L such that L > \t\ 7 ^ e > 1 one has for all i,j S L that 

\CiA<^, \ ai \ < \\(3\\ij07w- (29) 

This bound shows the exponential smallness of the matrix entries of C = e~ ltA in 
the distance d^j away from the main diagonal at a given time. Also, Eqs. d29l signify 
bounds needed to prove convergence due to mixing within the causal cone for large 
times. The proof will be given in the appendix. 
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4 Quantum central limit theorems 



Lindeberg central limit theorems are central limit theorems that do not rely on an iden- 
tically distributed assumption concerning the considered random variables. Here we 
formulate a quantum version thereof, which will in the main theorem be extended to a 
quantum version of a central limit theorem that in addition does not rely on independent 
random variables. We start by noting a number of useful facts on quantum characteris- 
tic function in order to proceed to formulate the quantum central limit theorem itself. 



4.1 Some facts about characteristic functions 



In this section, we will relate quantum characteristic functions to classical characteristic 
functions, by invoking Bochner's theorem. For a similar correspondence of quantum to 
classical characteristic functions, see, e.g., Ref. 11261 . For a fixed vector of phase space 
coordinates a <E <C L consider the function 6 : M — > <D, 



<f>(r) = ti[g Q D(ra)} = tr[£ c^ a >] = (c ri ^). 
One finds ^ (0) = for odd n G N and 

0( 2 >(O) = a(a), 0< 4 >(O) =/(ct), 



(30) 



(31) 



where (jr- n > denotes the n-fh derivative of <fi. Using the Cauchy-Schwarz inequality, one 
finds for even m £ N the following chain of inequalities. We suppress the argument a 
of b defined as in Eq. (fT3l ). use the fact that b< = —b and write go = ^ n £>n|'0n)(' ! / , n| 
for the spectral decomposition of the state g . We hence find for the m-th derivative 



\^ m \r)\ <Y,^\^\b m/2 e rh b m ' 2 \^ n 

n 
n 

<Y.Sn\^n\{^) m b m \^ n ) 



n 

1/2 



|0M(O)| 



(32) 



1/2 



^(3m)( ))l/2_ 



Furthermore, 0(0) = 1, and <fi is positive semi-definite and continuous at the origin 
(27). It follows from Bochner's theorem that is a classical characteristic function of 
a classical random variable. Hence, if (j)^> (0) exists, </> is twice continuously differen- 
tiable on the entire real line, i.e., as 0W(O) = 0, one has from Taylor's theorem 



— II < 



Hence, 



Ax\4> (2 \x)\{\-x) < 

<!>(«)>- ii < M?oi 



l^ (2) (0)| 



(33) 



(34) 
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for all a G <C L as a was arbitrary. 

If (/)( 4 )(0) exists, is four times continuously differentiable on the entire real line, 
i.e., one has from Taylor's theorem 



n)-i- 



(35) 



4.2 A Lindeberg-type quantum central limit theorem 

We now turn to the actual quantum instance of a Lindeberg-type central limit theorem. 
Let Ai, ... , A n C C be mutually disjoint sets. These take the role of independent but 
not necessarily identically distributed random variables. Consider 



Zi = 4>i{l) = (cxp[b(a Ai )}), i = l,...,n. 



(36) 



If|zj — 1] < 1/2 for all i, a complex logarithm is defined by the Mercator series. We 
have 



log(*) - ^ (2) (0)/2 < |log(zi) - (zt - 1)| + ^ - 1 - 4>T' (0)/2 



(2), 



< \\og{zi) - (zi - 



where 



(37) 



|log(zO - (z, - 1)| = |log(l + (* - 1)) - (* - 1)| < £ 



Ift -11" 



1 n + 2 ~ 1 1 ^ ?i + 2 



< \Zi - \\ z < 



l^ 2) (o)l 2 < l# } (0)| 



(38) 



Hence, 



( 4 )/ 



Now, let x,y,z € R. We set out to show that \e x+ly — e z \ < \x — z + iy\e" 
This holds for x = z, i.e., we may let x ^ z w.l.o.g. Then 



(39) 

<-{x,z} 



i z (e x - z+iy -l)=e z {x-z + iy) dtc t( - x ~ z+iy \ (40) 



i.e., using the mean value theorem, there is a c between x, z such that 

"i 



c z \ < c z \x- z + iy\ / dte t{x ~ z) = \x - z + 



(41) 



\x — z + iy\e c < \x — z + iy\e n 
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Now, we recall that \zi — 1| < 1/2 and \zi\ < 1, i.e., 

log 1^1 = -| !og|zi 
Hence, using the fact that = — \a^ ( |, we find 



(42) 



n 



< 



< 



]Tlogz. t -]T 

i i 

^log^-^] 



2 

q-^ti 

2 



3 max{- £\ | log |zi||,— J2i \"Ai 1/2} 



(43) 



which yields the following theorem. 

Theorem 1 (Lindeberg quantum central limit theorem) Let A4 C L, i = 1, . . . , n, 

be mutually disjoint sets, \<JAi \ ^ 1» and let < 00. Then 



JJ0(a A ))-e^ 



7 ™ 



(44) 



A dynamical quantum central limit for initially uncorrelated states qq = ®f =l Qi fol- 
lows directly from this theorem. The subsequent technicalities are due to the fact that 
we allow for correlations in the initial state. 



5 Main theorem 

We now turn to stating the main theorem. We arrive at a statement showing local 
convergence to a Gaussian state in trace-norm, under the weak assumptions on the 
initial state specified above. We first introduce a variant of a Bernstein-Spohn blocking 
argument, and then make use of the above quantum Lindeberg central limit theorem, to 
show convergence to a Gaussian characteristic function in phase space. We then show 
that this implies closeness of q$ (t) to a Gaussian state — so a maximum entropy state 
for given second moments. 

5.1 The Bernstein-Spohn blocking argument 

We divide — for a given time — the lattice into two parts: Into sites inside the causal 
cone, where mixing occurs, and sites outside the causal cone, the influence of which 
is exponentially suppressed. Inside the cone, we in turn split the region into several 
subsets, the role of which will soon become obvious. This choice of a blocking is 
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derived from the "room and corridor" argument on classical lattice systems presented 
in Ref. GJED. To this end let 



> b > 1 such that n := 



8e|i| 



> 1. 



(45) 



Then 



8e|t| > n(a + b) > 4e\t\. (46) 

Now define the set 

T = {j e L : dj, s > n(a + b)} (47) 
as the sites that are significantly outside the causal cone, compare Lemma [3] Define 
also the sets A = U"=i A an( l & ~ U™=i wnere 

A = {:) e £ : (i - l)(a + 6) < d JiS < i(a + 6) - 6} , 
ft = {j £ I : i(a + b) - b < d^ s < i(a + b)} . 



(48) 



Then £ = AU BU T — the entire lattice is covered by these sets — and S C A\ . What 



is more, An B = AnT = T n B = Ai n Aj 



B, n Bj 



for i ^ j. The 



sets A and B cover hence the complement of T in an alternating fashion, the first set 
Ai containing the subset S of interest for which we aim at showing relaxation. The 
goal is now to choose b sufficiently small (such that the influence of sites in B is small 
due to the large separation of the £>;), but large enough such that the Ai are sufficiently 
far apart (and therefore approximately independent) and the Lindeberg central limit 
theorem yields closeness of Xg s {t) to a Gaussian function in phase space. We will see 
that for sufficiently large |t|, this is achieved by the choice 



\t\ 2 / 3 
log |*| 



U|(2-Mi/(l+Mi))/6 j 



(49) 



for which one can prove the following Lemma. 

Lemma 4 (Correlations) Let Assutnptions\I]and\2\hold and let t, L be such that 



where ^ = fix/ (6(/ii - 
UP/ 3 



L 6/7 > \t\ > 2, 
1)). Then, 



log |*| < I* 



1/3+M 



a> b= It] 1 / 3 -" > 1, 



log |* 

and for A B, T as above, one has 

WtIWb,t\,Wa,t\ < c 2 \\(3\\i2- Mtl , 

k B |> AS |<16ec 2 |&S||!/3|' 2 ^ 



n > 1, 



(50) 



(51) 



CT£ 



1=1 
n 

E 

i=i 



< 



1 |t|M 

c2||/3||?(|p+80e|&S| 



log 1*1 



(52) 



|/Al<C4«||/3||?(J^ 



2\dS\ 
log |*| 
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The proof is again presented in the appendix. Here, the constant jii is defined as 
specified in the assumptions on the initial state, C2, C4 as in the Lemmata following the 
assumptions, and 

55 = {1 e S : d hC \ s = 1} (53) 

is the boundary of S. 



5.2 The main theorem 



Let A, B, T as above and write Da = D{a.j\) and similarly for other the sets. We 
have (we employ the triangle inequality three times) 



?c/2 



< 



(D C ) - (Da) 



n^)-ii e 



"aJ2 



(54) 



The first term is bounded by (we write again go = Qn\'4>n){4>n\ for the spectral 
decomposition of qq and employ the Cauchy-Schwarz inequality) 

(D A (D BUT - <J2en\(lpn\DA(D B UT-~L)Wn}\ 
n 

<Y,V^ (sn(^n\{D BuT - 1)^ (D BuT - iMn)) 



1/2 



< [{(D Bur -iy(D BUT -i)) 



1/2 



(55) 



m[((i - d but )) 

< (Wb\ + Wt\ + 2\(7 B , T 



1/2 



< Kur| 1/2 

1/2 



To bound the second term in Eq. ( f54b we let I = {m G N : m even, 1 < m < n}, 
J = {1, . . . , n}\I. Then, using the fact that d^^Aj > b for i ^ j, we find under 
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Assumption^ 



\\»A \[l>A 



< 



C5 

< — r- 2 h 

- frl/2 + ^2 



;n^n^)-(ii^)(n^ 

■iei jeJ iei jeJ 



;n^>(n^>-n<Ao 

iei jeJ i=i 



;n^><n^>-<n^>n^ 

iei jeJ iei jeJ 



(iP*>n<^o-n<AA« 

iei je.J 



i=i 



< 



C5 



51/2+M2 

and, using the fact that d,A t ,Aj •_ a f° r hj £ ^> i 7^ i> 



;IPa>-I1^> 



;6/ 



< 



< 



C5 



a 1 /2+^2 



;n^)-(^ 2 )( n 

ie/ 2^iei 



; n ^> - n <^o 
( n D Ai ) n <^o 

2/iei i^iei 

1/2+M2 



<-<m- 



This means that 

n n 



i=l 



< 



< 



c 5 n 



51/2+ M2 a l/2+M2 



C5 



8c 5 e 



(log|£|) 3 / 2 +^ 
|i|W3 



1^(1/12+^/6) 

which bounds the second term in Eq. d54l i. 

Now, for 1,1/ € 1 there is a i-q between a; and y such that 

e* - e v = [x - y)e x °, i.e., \e x - e v \ = \x - y\c x ° <\x- y\e max{x ' v} . 
Hence, for the third term in Eq. d54l ) we find (we recall that e>-£, a Ai < 0) 



< 



/2- 
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If \cAi I < 1> trie l ast term i n Eq. (EH may be bounded using the Lindeberg-type central 
limit theorem and we finally have 



< (\o B \ + \a r \ + 2|cr e , r | 



,1/2 



CT£ 



1=1 

C5 



(61) 



i=l 



^(1/12+^2/6) 



8c 5 e 



(log|t|) 3 / 2 +^ 
| t |a/te/s 



where bounds on the moments may be read off from Lemma |4] This yields the main 
theorem. 

Theorem 2 (Convergence of characteristic functions) For all C = {l,...,L} c 
N, all S C C, and every go fulfilling Assumptions\l§3\there exists a Fs(/3,t, L) > 
such that 



X§ s (t)(P) 



<F s (f3,t,L). 



(62) 



The function F$ : C' 5 ' x 1 x N -} B has the following properties. There is a 
recurrence time t rec > 0, which increases polynomially with L, and a relaxation time 
treiax such that Fg is decreasing in time for all \t\ £ [t re iax,trec\, and all S C {1,...,L} 
with 

T/'f/te Zaffer /zo/cfc a«<i |<S| rfoei «of depend on L, one has in particular 



lim lim F s (/3,t,L) = 0. 



(64) 



This theorem indeed shows convergence of characteristic functions of reduced 
states associated with some sublattice S to Gaussian characteristic functions in phase 
space. Note that the size of this sublattice may even slowly grow in time; one would 
still end up with a Gaussian characteristic function. Surely, in any finite system there 
will be recurrences, but with increasing L, these recurrence times i rec can be made 
arbitrarily large. 

The characteristic function 



corresponds to the Gaussian state qq if) with second moments 



h(t)kk = mMt)) + Mt)bl(t)) = [c-^ 7 (0)O, fe 



and one hafl 



tr[e G (t)Sj6i]=tr[&(t)Sj6i] 
[7(*)k< ~ 1 



< ||e-^ 7 (0)c ItA || = || 7 



(65) 



(66) 



(67) 



Here, 



| denotes the spectral norm. 
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Note that, if the initial correlations are translationally invariant, 7(0) and A commute, 
i.e., the second moments are conserved, 7(4) = 7(0). By virtue of Lemma[5] which 
we prove in the next section, trace-norm convergence of states is also inherited from 
this main theorem. We state this corollary for a sublattice S that is constant in time. 

Corollary 1 (Convergence of reduced states to Gaussian states) Let \S\ and ||7(0)|j 

be bounded by an absolute constant. Under the assumptions ofTheorem\2\andfor any 
e > there exists a recurrence time t rec > 0, which increases polynomially with L, 
and a relaxation time t re i ax such that 

Us(t) - QG(t)\\tr < e, (68) 
for all \t\ G [t r eiax, tree], where Qc(t) is the Gaussian state with characteristic function 

This means that in the sense of trace-norm closeness, local reduced states become 
Gaussian — and hence maximum entropy — states. 



6 Convergence of quantum states 

In this section, we will show that pointwise convergence of characteristic functions is 
inherited by trace-norm convergence on the level of quantum states. This argument 
generalizes the one of Ref. |27l — showing weak convergence for states from pointwise 
convergence of characteristic functions — and of Ref. [28] - relating trace-norm and 
weak convergence. The central idea is to consider two sets: On the one hand, this is 
the subset of phase space that supports, in a sense yet to be made precise, "most" of the 
characteristic function. On the other hand, this is the subspace of Hilbert space with the 
property that the projection of the state onto it is a positive operator that eventually has 
almost unit trace. One can then bound the error made when neglecting the complement 
of these sets and in turn relate the relevant sets in phase and state space. Inside these 
relevant sets, one can relate weak convergence to pointwise convergence of character- 
istic functions and use the equivalence of norms. Here we present an argument going 
beyond the results in Refs. l27ll28l and use a gentle measurement Lemma of Ref. l29l . 



6.1 Preliminaries 

To prepare the argument, consider for a single mode the upper bound to matrix elements 
of Weyl operators 



[n\D(a)\m) 



< e" 



M 2 / 2 uil"-™ 



min{m, n}\ 
max{m, n}\ 



|\ 1/2 



I n — rn I 



i{m,n} 



(M 2 ) 



(69) 



= :e"H / 2 p„,m(H 2 ), 



where a E (D and L„ are associated Laguerre polynomials, i.e., p n .m(\ct\ 2 ) is a poly- 
nomial of degree min{?7i, n} + \n — m\/2 = (n + m)/2 in \a\ 2 . Here {|n)} n= o i i j ... is 
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the standard number basis of the Hilbert space associated with a single mode. We can 
make use of this in the situation of having N modes: Define for some T > 

M := {a G <E N : \\af 2 > T} , (70) 

as the region far away in 2-norm from the origin of phase space and A4 := (C N \A4 as 
its complement. For a continuous function / : <D N — > <D with |/(a)| < 2 consider the 
bound (we write p n ,m(a) = ri l Pn«,m i (|a l | 2 ) and \n) = |m • ••n N )) 



daf{a)(n\D(a)\m) 



< da |/(a)| (n\D(a)\m) 

M 



da \f{a)\ (n\D(a)\m) 

M 



<2 da 
'm 



[n\D(a)\m) + \M\ max \f(a)\ 



< 2 [ dae-^/ 2 Pn . m {a) + \M\ max \f(a)\ 

Jm a£M 

< 2e- T / 4 / dae-IHI*/V,m(a) + \M\ max |/(a) 



-T/4 



+ \M\ max|/(a)|. 



aeM 



(71) 



These findings may be summarized as follows. For all N,i%i,mi £ N and every e > 
there is a cq > and a compact set C C (D w such that 



da f(a)(n\D(a)\m) 



< e + Co max |/(a)| 



(72) 



for all continuous / : € N -> C with \f{a)\ < 2. 
Now, for M C N w let 

P M =Y J \n){n\- (73) 

This will take the role of a projection onto the relevant part of state space. For states 
gi , Q2 one finds (see the Appendix for a proof) 



llei - &||tr < ||ex - PmQiPm + PmQiPm - &||tr 

+ |M| 3/2 max \(n\(gi - g 2 )\m) 

What is more, using the Heisenberg-Weyl-correspondence, one has 



In (gi - Q2)\m) cx 



da [Xei( a ) ~ Xq 2 (<*)] (n\D(a)\m) 



(74) 



(75) 



I.e., combining Eqs. (172174b . one has that for every N G N, e > and M C there 
is a constant cq > and a compact C C (D such that 



\§i - hWa < e + c max|xe 1 (a)-Xe 2 (o:)| 



\Qi - PmQ\Pm + PmQ^Pm ~ &||tr 



(76) 
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for all N— mode states gi, g 2 . 

Now, from the gentle measurement lemma of Ref. [29] it follows that 

\\g-PMgPAi\\tr<2(iv[g(l-P M )]y 2 (77) 

and for hi = b\bi one finds (we write M = N JV \M) 

N N N 

Nma,xtr[ghi\ > ^tr[^n,] = y~] rii(n\g\n) + X! X! n i( n \Q\ n ) 

i— 1 2—1 n£LM 2—1 n£M 

N 

> tr[f>(l - P M )} min V 7i, > tr[£(l - P M )]m, 

nEAI ^ 



(78) 



where the last inequality holds for the particular choice of 

AI = {ne N w : ^ n, < m}. (79) 

i 

Hence, to summarize, for every to, N £ N and e > there is a constant cq > and a 
compact C C <D W such that 

Mi - &||tr < e + comaxlxg^a) -X@ 2 (a)| 

[N\ l/2 f \ (80) 

+ 2 ( — J l^max Vte^A] + max vtr[^2^] J 

for all iV— mode states £ l5 g 2 - 

6.2 Trace-norm bounds 

We summarize the above findings in the following Lemma. It shows that closeness in 
trace-norm is inherited by pointwise closeness of characteristic functions. The last line 
corresponds essentially a bound to the local energy per mode, which can in all relevant 
settings be bounded from above. Note that all constants on the right hand side of the 
inequality can be made explicit. 

Lemma 5 (Trace-norm convergence from pointwise phase space convergence) For 

every e, Ci > and JVgN there is a constant cq > and a compact C C (D N such 
that 

Ui - faWtr < e + comaxIxgiCa) - X§ 2 (a)| (gn 
for all states g\, g 2 on a collection of N sites with, for all i, 

tr[gihi], tr[g 2 hi\ < c x . (82) 

Using Theorem |2] the pointwise distance between characteristic functions and e 
may hence be subsumed in e to yield the trace norm estimate in Corollary [TJ 
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7 Summary, discussion and relationship to notions of 
typicality 

The above theorem shows that weakly clustering states locally relax to maximum en- 
tropy states. This holds true in the strong sense of trace-norm convergence, to states that 
have maximum entropy for given second moments. Physically speaking, this means 
that locally, the system appears as if it had entirely relaxed, without the need of an ex- 
ternal environment. Every part of the system forms the environment of the other, and 
merely mediated by local interactions, the system relaxes. This is a proven instance 
of a situation in which under non-equilibrium quenched dynamics, closed systems do 
locally relax. As has been mentioned before, it can hence be viewed as a proof of the 
"local relaxation conjecture", generalizing the findings of Ref. J7), and constituting to 
the knowledge of the authors a first proof of local relaxation in a situation not requiring 
product or Gaussian, meaning quasi-free, initial states. That is to say, the statement 
shows that for a large class of initial states, locally, the system quickly looks relaxed 
and stays like that for long times. 

We finally mention a potential link of the above results to the concept of "typi- 
cality": The term typicality refers — roughly speaking — to the idea that when quan- 
tum states are randomly chosen from a suitable ensemble of large-dimensional pure 
states, almost all realizations will share properties like reduced states of subsystems 
having large entropy, or being close to a Gibbs state. This effect is a manifestation of 
the phenomenon of concentration of measure. The simplest such situation is realized 
when drawing random pure state according to the Haar measure on T~L — (D <g> <D D : 
Then the von-Neumann entropy of the rf-dimensional subsystem will, for large D, be 
almost certainly be almost maximum, with any deviations being exponentially sup- 
pressed in probability. Similarly, if states are drawn from an ensemble of pure states 
that are eigenvectors of a Hamiltonian with eigenvalues close to a limiting point in the 
spectrum, defining an energy, then one can specify technical conditions under which 
one can prove that almost certainly, subsystems will locally look like the reduction of 
a Gibbs state to that subsystem |fT71 [TBI . That is to say, almost all states appear to 
have high local entropy. This is strictly proven only in case of finite-dimensional con- 
stituents, although it is to be expected that the same intuition carries over also to the 
situation of infinite-dimensional constituents. In the light of this observation, it appears 
plausible that starting from an "untypical initial situation", systems will eventually be 
driven to a "typical one" by virtue of unitary time evolution. 

Needless to say, time evolution defines merely a one-dimensional manifold, pa- 
rameterized by time, and from the above argument alone, there is no guarantee that 
the system will indeed arrive at a typical situation (although steps can be taken that 
show that for most times, albeit unknown ones, the system appears locally typical |[T9l 
in a non-equilibrium situation). Here — for a class of models — a significantly stronger 
statement is proven, in that after a short time, the system will certainly come arbitrar- 
ily close to a maximum entropy state under the constraints of second moments and 
will stay there. It would be interesting to take steps to combine such kinematical and 
dynamical approaches towards understanding equilibration in closed quantum systems. 
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9 Appendix: Proofs 

9.1 Lieb-Robinson bound (Lemma© 

Under periodic boundary conditions, the entries of C take the form Cij = Ci-j — 
Cdi } , 

1 L 

Ci = - V e ^iki/L e -itx k ^ Afc = _ 2 C0S (2irk/L). (83) 
Li ' 

k=l 

As Ai j = for dij > 1, we find from Taylor's theorem the estimate 



fe=l 



< max 

feG£ 



E 



(-itXk) r 



\t\k\ d ^ , \2t\ d ^ 
< max ' < < 

lis £ £ Cby J ■ J ' 

i.e., for A C C we find from Eq. ( fTTT ) that 



2e[t| 
4j 



(84) 



En < EE^-ii^i ^ li^iii^E 

ieA ieAjeS J i£A 



2e\t\ 



< ||/8|| i max ^ 



|/3||i max > 
1 11 jes ^ 



2e|t| 

cU,S 



/2e|t[ 

V°U<S 



1> 



ie.4 



< 2||/3|| 



2e|t| 

dA,S 



dA,S oo 



E 

d=0 



( 2e\t\ 

\ d A,S 



<4||/3||i2 



-<2.a,s 



(85) 



where the last inequality holds for 4c|t| < dj\,s- 
Furthermore, 

L L 

Ci = ±J2e 2Hcos{2 * k/L) e 27rikl/L =: -^.g(27rfc/L)c 2 ^/ i ! 



k=l 



k=l 



(86) 
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where g(cj)) = G 2ltcos W can be written as 

oo 

3(0)= E 9ne in t 



n— — co 



<,)n 



1 f 

d0 ff(0)e -in * = — / d<£e 2itcos ^e i7 ^ = iV„(2t) (87) 



27T J 27T Jo 



2iT 



d0e in *Q2e 2itcosW , 



"2^F7 ^ ^ 
where J n is a Bessel function of the first kind. Hence, for \t\ > 1 

\ 9n \ < A r^\ d y—^\ < M±^l < 6 M! (88) 



27m 2 J n 2 n 2 

and 



oo L oo 

3 = 1 E 3nE e2Tifc( " + ° /L - E fo^eLi 

n— — oo fc— 1 n— — oo 

r2ir 



z=l z=l 



i.e., for L > |i| 7 / 6 > 1 (we use a bound obtained in Ref. Il30l ). 



C dl \ < 



1 /' 27r 



2^ Jo 



(2*)| 



< 6(^-4)1^ _^ l < 37 

l 2 i" \2t\m - jipTs- 



9.2 Moments 

We will need the following bound. Let e,fi > and A C C. Then 

1 °° 1 

£[l + ^] 2+ " +£ ^ [l + rr[l + r] 1+ "S' iM ' 

2C(14 



(89) 



= E 0J«-< = a" / d( M0)e" + E + fli-+0 . 

Z = -00 "' 2 = 1 

where, for \t\ > 1 and < I < L/2, which we can assume as Cij — C^ ., 

oo oo 

E + .^1 < * E (l^W + = (90) 



(91) 



(92) 
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9.2.1 Second moments (Lemma [T) 

In the following we assume that t and L are such that L > \t\ 7 / % > 1. Under Assump- 
tion Q] we find from Eqs. d 1 8192b and Lemma [3] that 



^B|<El°* a J-|(|^ S i)| + |M) 



< 



ieA 

jeB 
37 2 



|i| 2 / 3 



|/3||?^(2 |<StS J ->| 



ieA 
jeB 



< 



37 2 (2 Cl + l) 2 



1 



(93) 



| t |2/3 ll/ 3 H?E [1 + rfij] 2 + M 1+ei 

jeB 

< 37 2 (4 C1 + 2) C (1 + eJMl m " {lAm L 

|t| 2 / 3 [1 + <Ub] 

for all i,Bc£. Similarly, for all A C C, i = 1, . . . , n, with A l n A, = for i ^ j, 
we find (.4 := \J. Ad 



Ei^i<^^ii/3ii?EE E [1+d J^ 



=1 fce^ ze-AVA; 



(94) 



<c 2 ||/3|| 2 £ 



I Al 

\t\*' 3 [l+d Ai , A \ Ai 



Now, for all A, B C £ with .4. n $ = one has cr^us = + cfb + 2(7.4, i-e-, 

n n 

f£ = f£V4 + 2&£,\A,A + < J A = °C\A + %Oc\A.A + E '- 4 ' + E <7 - 4 *.^' (95) 



i.e., 



< k£ul+2|a AA ^|+ C2 ||/3|| 2 ^ 



|t| 2/3 [l+<U\, 



(96) 



Finally, under Assumption Q] and for i, .4 such that 4e < 4e|t| < ci^s, we find 
from Eqs. ( 118192b and Lemma [3] that 



Wa,b\ < E 



(2 Cl + 1)1^^1 l __^37(2ci + l)||)3||i 

X + 



- > IcKil max > 



jeB 



ieA 



(97) 



< 37(4 Cl + 2)11/311^(1 + ex) £ |oj| < c 2 ||/3|| 2 2-^' 

ie.4 
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9.2.2 Fourth moments (Lemma |2) 

From Eq. (Qj]) we find 

\Ia\ < £ i^^«fe«H(l(&l6]^^>i + K6l&fe&J^)l + K&!^^)i 

i ,j,k,l<=A 

+ KSfcSJStS,)! + KS*S5S,SJ>| + K& fc $$&}>|) (98) 

<3 2 l« I « J «fea/|(2|(SlS]S fe S0|+44, J |(Si6 / )|+<5 fcj 4 z ), 
where, under AssumptionQ] we have 

^ (4c 1 + l) (99) 

([l+dr, s ]tl+<i t ,„]) 1+ei ' 

(r,s,t,«)eP(i,i,fc,0 

i.e., under Assumptions[T|and|2]we have 

IJUI < I<*«j«*«iI E ([i+^.^S+^'i-'^ - (100) 

ij',fc,ZG.A (r,s,t,u)£P(iJ,k,l) 

Now let (r,s,t,u) be a given permutation of (i,j,k,l). Then, using the geometric 
mean inequality and Eq. 



|ai«j||qfcQn| _ larjjasjjgtjlgg] 

([l+d r ., s ][l+d t ,„]) 1 + min < E i' e 2> " ([l+d 7 ., s ][l+<i t ,„]) 1+min{e i' e 2> 
i,j,k,l£A r,s,t,u£A 



L [l + dr, 3 ] 1 + minlei ' e2} 



r.seA 



2 



(101) 

< r v k-i 2 , y 

<4C 2 (l + min{e 1 ,e 2 })[E|a l 



|2 



i.e., using Lemma[3]and the fact that J2i&c \ ai \ 2 = ll/^lli an d \P\ — 24, we find for 

L > \t\ 7 / 6 > 1 



n n 

E l/A I < 96(6c 3 + 3(4ci + 1))C 2 (1 + min{ ei , e 2 }) £ [ £ 



i=i *ieA 



< 96(6c 3 + 3(4cj + 1))C 2 (1 + min{ei, e 2 })||/3|| 2 max V 



iax 2^ | ail 

jeA 

zV 2 m j. ^„\Mi/qii 2 ii/qii 2 max - 1- 4 -! 



< 96(6c 3 + 3(4d + 1))37V(1 + min{ ei , e2 })||/3||^||/3||fH^ 

(102) 
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9.3 Blocking argument (Lemma 31) 

In the following we write r = 8e|t|. 

9.3.1 Distances 

For all k £ A4, 1 £ Aj we have (we pick s £ S such that d; s = 5) 

(i-l)(a + b) < d k . s < d k . s < d k ,i + di, s = d k ,i+di, s < d k ,i+j(a + b)-b (103) 
and (we pick s £ S such that dk. s = dfc.s) 

(j-l){a + b) < di, s < d t>s < di,fc + d fc , a = d Lk + d k , s < d k ,i + i(a + b) - b, (104) 
i.e., 

dAi.A, = min d k l >\i — j\(a + b) — a. (105) 

keAi 

Furthermore 

T 

ds,T = yiinds,; > n(a + b) > — . (106) 

9.3.2 Cardinalities 

For i £ C\S pick s, £ dS such that di t s = e?i. Si = di : Qg- Then 

\M = E E = i 5 i+ E E 6 <s.i = \s\ + E E 

0<l<ai£C l<l<aieC\S l<l<ai£C\S 

^\ S \ + E E E **..,i<|S|+2|aS|a. 

l<i<asG9SiG£\5 

(107) 



Similarly, for i ^ 1 



\Ai\<2\dS\ £ l<2\dS\a (1Q8) 

(i-l)(a+b)<l<i(a+b)-b 

i.e., 

n 

\A\ = \A\ < \S\ + 2\dS\na < \S\ + 2\8S\t. (109) 

Furthermore, 

n , 

|£|<2|&S|^ J2 !<2\dS\nb<2\dS\ — . ( no) 

i=l j(a+h)-6<Z<i(a+b) 
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9.3.3 Moments 

In the following let t and L be such that 

L 6/7 >\t\>2, \og\t\ < \t\^ 3+fi , 
where /j, = /Lti/(6(/Xi + 1)). Then 



|f[ 2 / 3 
log |t| 



= a > b= \t\^ 3 -^ > 1, 



8e|t| 



= n > 1. 



(Ill) 



(112) 



Combining the above bounds on distances and cardinalities with Lemmas Q~|[2] we find 
under Assumptions 1 1 12l that 



Wt\,Wb,t\,Wa,t\<02\\P\\P- 4 * W , 



and 



M,|cr AS |<2c 2 |aS|||/3||? 



Tb 



16ec 2 \dS\\\f3\\i 



2 |t|V3 6 



a|t| 2 /3 

= 16ec 2 |S5|||^||?|i|-^log|t|, 

< laeurl + 2|<7 Bur ^| + c 2 ||/3||? ^p^n 

< ksl + krl + 2|tr Sl r| + 2|tr B ^| + 2|ct at | 
+ c 2 ||/3||?(J^3+16e|S.S|i^ 

< 5c 2 1|/3|| ?2- 4 °l*l 

r 3|t|^log| 

c 2 ||/3||?(^L + 80e|a5||t|- /1 log|i|' 
El/Al<c 4 ||/3||^||/3||?^< C4 ||/3||^||/3||?^ 



■c 2 ||/3||f ( + 16e|55| 



< 



f'4 



|*| 2 / 3 

ll/Wlll 



| * | 2 / 3 



2\dS\ 

tp/3 "i" EiT?T 



9.4 Closeness of quantum states 

We have for any states Q\,Q2, 



\QI - &||tr <\\h - PmQiPm + PmQ-zPm - (02 1 1 tr 

+ \\PmQiPm - Pm&PmWk, 



(113) 



|t|V3 



(114) 



(115) 
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treating the relevant part of state space. Within that relevant part, we can make use of 

\\PmQiPm - PmQ2Pm\\u = \\Pm(qi - h)PM\\tx 

= tr 



1/2' 



(Pm(qi - Qi)Pm{qi - h)PM^j 
~y]{n\ (Pm(qi - Q2)Pm{qi - Q2)Pm} \n 



1/2 

(116) 



giving rise to the bound 



\ 1/2 

\PmQ\Pm - PmQ2Pm\\\x < 2J ( { n \ p M{Qi - Q2)Pm{qi ~ (h)PN\n) 



n \ m / 



1/2 



(117) 



£ ]T |<n|(&-&)|m; 



\ 1/2 
|2 \ 



nGM VmeM 



< |M| 3/2 max |(n|(gi - &)|m>| . 
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